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Abstract 

The known ratio of the branching fractions for D* — > D s tv° and D* — > D s y may be used to extract 
the coupling of 77 mesons to strange quarks once the value of the n — n mixing angle is known. This 
requires that realistic models for the spectra as well as the magnetic dipole (Ml) decays of the heavy- 
light (Qq) mesons are available. The coupling of n mesons to light quarks may then be estimated 
using SU(3) flavor symmetry. Applied to the quark model for the baryons, an nNN pseudovector 
coupling constant of f V NN = 0.35 1q 25 ^ s obtained. If the charm quark couples significantly to the 
n meson, as is suggested by the decay mode ip' — > J/ip n, then somewhat larger values of f n NN can 
be obtained. These values are sufficiently small to be consistent with phenomenological analysis of 
photoproduction of the n on the nucleon and the reaction pp — > ppn. 
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1 Introduction 



The flavor symmetry breaking D* — > D s ir° decay mode is mainly due to a small isoscalar 77 meson 
component in the physical ir° meson as only the 77 meson can couple to the strange quark in the strange- 
charm mesons [Q, Q]. As the magnitude of the 7r°— 77 mixing angle has been at the focus of a considerable 
amount of theoretical and phenomenological analysis, its value is now rather well known B]. The 
expression given by chiral perturbation theory to lowest order H is, for small m , 

m = -. " " - 0.010 (1) 

4 m s — mo 

where tuq is the average of the u and d current quark masses. To this should be added the contribution 
from ir° — r)' mixing, which is smaller because of the larger mass of the 77', and which may increase the 
value of 9 m by at most about 30% || |(| . The resulting effective n — 77 mixing angle is therefore about 
- 0.012. 

Combination of the notion of the spontaneously broken approximate chiral symmetry of QCD with 
SU(3) flavor symmetry implies that the coupling of the baryons to the octet of light pseudoscalar mesons 
(ir, K, 77) takes the form 

Cnnm = i-^f- 4>n Isl^dn m a X a ip N , (2) 

where ^n denotes the baryon octet field, m a the meson octet field, A a the SU(3) Gell-Mann matrices, 
f m the decay constants of the mesons, and gA the axial coupling constant of the baryons. The current 
experimental value of gA as determined from low energy nuclear beta decays is gA = 1.2573 0, while 
the phenomenological value is gA = 1.267 ± 0.004 ||. If these values are employed together with the 
Goldberger-Treiman relations for nucleons, then "effective" pseudovector coupling constants may be 
defined as 

m n 1 m„ 

Unn = 7rr9A^ 0.95, frjNN = gA ^ 1.79. (3) 

V O ^Jri 

The above value for the irNN coupling constant agrees completely with that used in the most recent 
realistic phenomenological nucleon-nucleon interaction models, the parameters of which have been de- 
termined by fits to nucleon-nucleon scattering data ^TJ]. Such fits do not, however, constrain the 
value of frjNN very well, and therefore its value remains uncertain. Different realistic phenomenological 
potential models give values between ^0.5 and ^3.6 for this coupling [ pT| , while analysis of data on 
photoproduction of the 77 on the nucleon indicate that the value for J^nn should not exceed 0.64 p"2[ , 
a value which also seems to be consistent with data on 77 production near threshold in proton-proton 
collisions [jl3| . 

The point of the present paper is to show that it is possible, given the value of the 7r°— 77 mixing angle, 
to estimate the strength of the coupling of the 77 meson to strange quarks from the empirically known 
ratio of the branching fractions for D* — * D s ir° and D* — > D s j. This requires that realistic models 
for the respective hadronic matrix elements and Qq wavefunctions are at hand. From the 77SS coupling 
obtained, the corresponding r/qq coupling may be estimated using flavor SU (3) symmetry. One may 
then, by standard quark model relations, obtain useful and constraining information on the 77— baryon, 
and in particular, the 77— nucleon coupling f n NN- The main source of error in this analysis is the rather 
large uncertainty in the empirical determination of the 7r° and 7 branching fractions for the D* meson. 
This error is likely to be much larger than that introduced by SU(3) breaking effects into the current 
analysis. 
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If one considers the coupling of the octet of light pseudoscalar mesons to the light (u, d, s) quarks, 
one obtains analogously to eq. (0) the couplings 



A 



£qqM = i-^T- ^9 75 7 A m aK Ipq, 



(4) 



where g q A denotes the axial coupling constant for constituent quarks. For the pions and the rj meson, 
the empirical values of the decay constants are f v = 93 MeV and f n = 112 MeV, respectively, so at 
least in this case the SU(3) flavor symmetry is broken only at the 10% level. Combination of the chiral 
coupling (W) with the representation 



m a Xa = V2 



V2 



V 



7T 

K~ 



V2^ 



n 



K+ 
K° 
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d 



leads to the following definitions for the quark-level pseudovector coupling constants f m qq, 
analogous to those given by eq. (||): 



(5) 



which are 



r q 

jTTqq — ~^J-gA' 



_ m v g 
Imq ~2^,f, 1 9A ' hs 



(6) 



The above relations then suggest that the magnitude of the coupling of r) mesons to u, d quarks should 
be one-half that of the rj coupling to strange quarks, independently of the r/ meson mass. The numerical 



value of g A is not very well known, although typical values fall in the range 0.87 - 1 [14, 
of the pionic decay widths of the D* mesons show that value range to be realistic [16 



m 

17| 



Calculations 
In the static 



quark model the meson-quark coupling constants eq. 
constants given by eq. (J3j> as (see e.g. ref. JlTy) 



(||) are related to the mcson-nucleon coupling 



NN 



' firqq, frjNN — fr) 



(7) 



Thus for gA ~ 1-26 one would predict a value of g q A ~ 0.75, which may be considered as a lower bound 
for the constituent quark axial coupling [ fl7| . If the relation (||) is employed with values of g q A in the 
range 0.87 - 1.0 [[u], [l^], then the corresponding values of f vss would fall between -2.5 and -2.8. Note 
that the value of g q A is conventionally extracted from the axial (spatial) current coupling, as the value of 
the effective axial charge coupling of confined constituent quarks depends on the form of the confining 
interaction Jl6| . 

The empirical observation of the charmonium decay modes ip' — > J/ip ir° and tp' ~ * J /4> r l indicates 
that charm quarks also couple to the ir°, and may potentially influence the D* — * D s ir° decay. The 
fairly large empirical width for ip' — * J/tpr] of ~ 7.5 keV which corresponds to a branching fraction 
of ~ 3% suggests that besides qq and ss the rj meson may have a cc admixture. While this would 
allow direct coupling of the rj meson to charm quarks, the charm components of the rj and rf mesons 
are phenomenologically, as well as theoretically, constrained to be rather small || [l8| [lj|. The known 
decay width for tp' — > J/tprj |^(| 21 and the Hamiltonian model for the cc spectrum of ref. j22| may 
nevertheless be used to explore the consequences of a nonzero 77CC coupling for the decay D* — > D s ir° 
and thereby the resulting value for J^nn- For this purpose the effective rjcc coupling is here taken to 
have the form (0). 
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This paper is divided into 5 sections. In section 2, the method for extracting the r]ss coupling from the 
ratio of the empirically known branching fractions for the decays D* — > D s tt° and D* — ► D s ^ is described, 
along with the models for the 7r° and 7 decays of the D* . It is shown that a relativistic treatment of both 
the 7 and ir° decays of the D* is called for. In section 3 it is shown that the two-quark exchange current 
contributions to the 7 decay amplitude, which are associated with the scalar confining and vector one- 
gluon exchange (OGE) components of the Qq interaction, provide a way of obtaining a realistic value 
for the matrix element for 7 decay, which is not possible in the impulse approximation. These two- 
quark contributions appear here because of the explicit elimination of negative energy components in the 
Blankenbecler-Sugar (BSLT) reduction of the Bethe-Salpeter equation. Section [| contains the numerical 
results for the ir° and 7 matrix elements, and an estimated upper bound for the ijcc coupling. Section 
5 contains a concluding discussion of the results obtained for the coupling of rj mesons to quarks and 
nucleons. 

2 The D* s -> L> s 7r° and D* s D s j Decays 
2.1 Extraction of the tjss coupling 

Application of the relations (|J) together with eq. (|J) yields a coupling of rj mesons to strange quarks in 
terms of fnss- The coupling of the ir° meson to the strange quark may then be expressed in terms of the 
"effective" mixing angle 6 m as , 

£ssn° = i-^Om Tpslbl^d^a lp s , (8) 

where sin(0m) has been replaced by 9 m , which is a very good approximation for values of 9 m ~ 0.015 
rad. In eq. (|8j), 6 m now corresponds to the sum of the ir°— r\ and ir°— rj' contributions as discussed in the 
previous section. 

The differential decay width of an excited heavy-light (Qq) meson may be expressed in the form 

dT = J_Ml lT ,2 (Q) 
dn 8tt 2 Mi 11 M ' 19 ' 

where q denotes the momentum of the emitted particle, T/j is the hadronic decay amplitude for the 
process in question, and Mt/Mi is a normalization factor for the quarkonium states analogous to that 
employed in rcf. jL7). The width for the process D* — * D s ir° may then be obtained directly from the 
expression for D^* — > D^ir given in ref. Jl6[ | by the replacement g q A /2J v — > /j7ss^m/w^, giving 

T{D* S - D s n") = \M„f, (10) 

D7r M d* rn^ 

if the 7T° emission takes place at the strange quark. Here M.^ is a radial matrix element for pion emission. 
Similarly, using eq. (^) the decay width for the radiative Ml transition D* — ► D s j may be written in the 
form |2^| 

16 M Ds 

where a is the fine structure constant. In analogy with eq. ([To|), M 1 is a radial matrix element for Ml 
decay, which includes the inverse quark mass factor (see below). By means of eqs. @ and Q, the 
ratio of the ir° and 7 decay widths of the D* meson is then obtained as 



TiD^D^^-^aq^M^ 2 , (11) 



(\M 7T I 



(12) 
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Note that in the above equation, the dimension of is [MeV] -1 . Through use of the empirical ratio 

of pion and photon momenta [^0) known to be approximately 139/48 and the 77 meson mass of 547 MeV 
one may solve for the coupling constant f nss to get 

As the ratio of the 7r° and 7 decay rates is experimentally known, albeit with quite large errors, to be 
0.062 ± 0.028 ^0|, it is, given the rather well known value of 9 m , possible to obtain an estimate for the 
coupling constant f 1)ss - This requires a model for the matrix elements in eq. (|l3]), as discussed below. 

Eq. (|l^) was obtained by assuming that the ir° decay takes place only at the strange quark. However, 
the observed ip' — > J/iprj decay shows that the 77 meson also couples to charm quarks, which implies 
that 7T° emission by the charm quark may also take place in the D* meson in addition to 7r° emission by 
the strange quark through 7r° — 77 mixing. If the pion emission at the charm quark is taken to proceed 
through a coupling of the type M), then the effect of tt° emission by the charm quark may be taken into 
account by replacement of eq. ( |13| ) with the expression 

R = C ^|VW 7 1 2 - 4.814 fm- 2 , (14) 

where R is defined as 

R=\f vss M% + f. ncc Ml\ 2 . (15) 

Here denotes the radial matrix element for pion emission by the charm quark. Because of the small 
momentum of the emitted 7r° and the heaviness of the charm quark, this matrix element will be very 
close to 1. The phenomenological implications of a nonzero contribution to ir a decay by the charm quark 
will be explored in section |5|. 

The coupling constant f vcc can only be estimated from the empirically observed decay ip' —> J/ipr] 
which, although somewhat suppressed by the orthogonality of the ip' and J /ip wavef unctions, is known pp| 
to have a rather large width of ~ 7.5 keV corresponding to a branching fraction of about 3%. However, 
reliable estimation of the coupling constant /, )cc from the decay ip' — * J/ip rj is difficult for several reasons. 
The first is that the matrix element for ip' — > J /ip 77 is very small in the nonrelativistic approximation 
because of the orthogonality of the wavefunctions. In a rclativistic calculation this matrix clement is 
typically increased by an order of magnitude. Furthermore, the axial charge component from a coupling 
of the type (||), along with the associated two-quark effects, may also contribute significantly to such a 
decay [jl6). Thus if the axial current component of eq. (||) is used with the replacement f VS s&m — > fijcc, 
then the value for f vcc obtained by fitting the empirical width for 77 decay of the tp' should be expected 
to represent an upper bound only. The expression for the width so obtained is 

w-w^w^r 9 " 1 ^ 12 (16) 

Here M v is the matrix element for 77 decay, which is given in the next section together with the matrix 
element M.-^. Note that in eq. (|l6|), the difference of a factor 8 as compared with eq. ( |l0|) arises from 
the product of a factor 2 from the spin sum for a triplet-triplet transition and a factor 4 from the sum 
of the quark and antiquark contributions to 77 decay. In section ^| it is shown that evaluation of eq. (16) 
leads to an upper bound for /^ cc of about ~ 0.8. 
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2.2 Matrix elements for pseudoscalar emission 

When the strange and charm constituent quarks are treated non-relativistically, the radial matrix element 
for 7T° decay, M.„, is 



M™= I dru a (r)jo(^)^ 1-^7^+0(0, (17) 



OC 



I 2 i" 24 



o 

where u(r) is the reduced radial wavefunction of cs system. Because of the smallness of (48 MeV) 
and (r 2 ) (~ 0.2 fm 2 ), the numerical value of the matrix element A4^ R is very close to unity. If the 
constituent strange and charm quarks are treated as Dirac particles, then the matrix element (0) for 
7T° emission by the strange constituent quark is modified to 

-i />oo />oo />oo />1 

M* cl = - dr'r'u(r') / dr r u{r) I dP P 2 dz '(P,q v , z) 
77 Jo Jo Jo J-l 



where the function F^, which describes the strange quark contribution to pion emission by the D*, which 
is a l 3 Si state, to the D 8 ground state from the Lagrangian (|J) may be expressed as jl6| 



(E> + m s )(E + m s ) ( P 2 gj/A 

t a [P,q^z)-\j 4EE , i{ E' + ms ){E + m ; ^- 



Here the energy factors of the strange quark are defined as E = ^/m 2 + P 2 + g 2 /4 — Pq^z and E' = 
\Jm 2 s + P 2 + q 2 /4 + Pq-nZ respectively. In the above expressions, P is defined as P — (p' +p)/2 in 
terms of the relative momenta in the center-of-momentum system. In realistic models for the heavy- light 
mesons jfj|, and the baryons typical values of the constituent masses for the strange and 

charm quarks fall in the range m s — 400 MeV - 600 MeV and m c = 1.3 - 1.6 GeV, respectively. Because 
of the small constituent mass of the strange quark the value of the expression ([l9|) is expected to deviate 
significantly from the nonrelativistic limit (|l7|). 

Because of the small value of the momentum of the 7r°, one may to a very good approximation, as in 
eq. dl7|), set q v to zero in the expression (|l9|), which then reduces to 

1 2 m 

F? (P) = Urn F? (P, q w ,z) = - + - . (20) 
9^0 6 6 hi 



In this approximation the relativistic matrix element for tt° decay may then be expressed as 

2 roo />oo />oo 

M™ = - dr'r'u(r') drru{r) dP P 2 F* (P) j (r 1 P) j {rP) . (21) 



Eq. ( |20| ) thus indicates that in the ultrarelativistic limit, the matrix element for ir° decay assumes the 
value 1/3. For m s = 560 MeV, the value of the matrix element ( pl| ) is ~ 0.8, which shows that the strange 
constituent quark has to be treated relativistically. In addition, two-quark mechanisms associated with 
intermediate negative energy quarks were shown in ref. [ fl6| to give a small contribution to the matrix 
element for pion decay. These effects will be discussed in more detail in section ^. 
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In order to obtain the matrix element for pion emission by the charm quark in the D* meson, it is 
sufficient to make the substitution m s —> m c in all equations of this subsection. On the other hand, the 
matrix element for the decay ip 1 — > J/ip rj takes the form 



M 



Rcl 



dr' r'uj/,p(r') I drru^>(r) 
in Jo 




/ drru r (r) / dP P 2 I dz F^{P,q v ,z) 
Jo Jo J-l 



16 



Pq v z 



(22) 



where Uji^ and u^> denote the reduced radial wavefunctions for the J/ip and ip' states, respectively. The 
factor F£ may be obtained by substitution of m s — + m c and — > q n in cq. (|l9|). 

2.3 Matrix element for 7 decay 



In the nonrelativistic approximation the matrix element for radiative Ml decay of a cs meson in eq. (13) 
is f§ 

1 r 2 1 

12 m c m s 



M 



NR. 
7 



(23) 



Although the photon momentum g 7 is somewhat larger (139 MeV) than q^, it is still small enough so 
that the Ml approximation may be considered valid. If the strange and charm constituent quarks are 
treated as Dirac particles, then the canonical boosts lead to relativistic modifications of the spin-flip 
magnetic moment operators |2j| . The relativistic matrix element for Ml decay may thus, in analogy 
with eqs. ( p"8| ) and (pl|), be written as 



M 



Rcl 



p CO 




I dr' r'u(r') 


/ drru(r) 


10 


Jo Jo 



dPP 2 F2(P)j (r'P)j (rP), 



(24) 



where the factor F? is defined as 



1 

12 



-P C {P) 



-fl(P) 



(25) 



The functions f?{P) represent correction factors to the static spin-flip magnetic moment operators, and 
may be expressed as 



3K. 



(26) 



for the charm and strange constituent quarks. Here the energy factors are defined as Ei = \J P 2 + mf. 
From eq. (23) it is evident that there is destructive interference between the contributions from the heavy 
and light quark currents to the Ml decay rate of a heavy-light meson. Insertion of standard values of 
m c ~ 1500 MeV and m s ~ 500 MeV for the constituent quark masses reveals that the light constituent 
quark contribution is somewhat larger than that of the heavy quark, which leads to an overall negative 
value for the 7 decay matrix element, cf. Table ||. However, if the relativistic form, eq. (|2~i|), is employed 
then the charm and strange quark contributions very nearly cancel each other. This occurs because of the 
larger relativistic suppression of the strange quark contribution. The end result is, that in the relativistic 
impulse approximation the width for Ml decay of a charged heavy-light meson is very small fl23f. 
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This result is, however, unrealistic, since two-quark mechanisms associated with intermediate negative 
energy quarks will contribute significantly to the matrix elements for Ml decay. These mechanisms, which 
are illustrated by the Feynman diagrams in Fig. [I] appear because of the explicit elimination of negative 
energy components in the Blankenbecler-Sugar reduction of the Bethe-Salpeter equation. The resulting 
two-quark magnetic moment operators have been calculated for equal constituent masses in ref . p7| . It 
has been shown in ref. p9[ | that the two-quark contribution associated with the scalar confining interaction 
is required for obtaining agreement with experiment for the Ml decays in the charmonium (cc) system. 
Without that contribution, the width for tp' — > J/ip 7 would be overpredicted by a factor ~ 3. 

However, the situation is much more complicated for the heavy-light mesons because of the uncertain 
structure of Qq interaction. Relativistic effects are also likely to be substantial for the two-quark operators 
because of the low masses of the light constituent quarks. It will be shown in the next two sections 
that once the two-quark effects associated with a Qq interaction formed of scalar confining and OGE 
components are taken into account, then the matrix elements for Ml decay of the _D ± * and Df* may be 
restored to the same order of magnitude as the non-relativistic estimate in eq. (|23|). The form of these 
two-quark operators will be outlined in section ||, and the numerical results are given in section |[ 



3 Two-quark contributions to D* decay 
3.1 Two-quark contributions to 7 decay 

In the reduction of the Bethe-Salpeter equation for the sc system to a Blankenbecler-Sugar equation, the 
negative energy components are eliminated in the impulse approximation. These do however contribute 
to the current matrix elements of the sc system as transition matrix elements, as illustrated in Fig. |} 
In the Blankenbecler-Sugar equation framework these transition matrix elements have to be included as 
explicit two-quark current operators ifjcfl . These two-quark matrix elements have been shown to give a 
large correction to the spin-flip matrix element of the single quark current [p9| . 




Figure 1: Negative energy Born diagrams for photon emission by a cs meson. The diagrams shown 
correspond to both time orderings of the 7 emission from the charm quark. Two similar diagrams 
describe 7 emission by the strange antiquark. The quark momenta are defined according to p a = p\ + q 
and pb = p[ — q, while &2 is the momentum transfered to the strange antiquark. Note that only the 
negative energy component of the intermediate quark propagator is to be retained. 
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The qualitative effect of this two-quark current may be understood in terms of a shift of the constituent 
quark mass by the scalar confining interaction m — > m + V c (r). Since the constituent mass appears 
in the denominator of the magnetic moment operator, this mass shift will lead to a reduction of the 
corresponding single quark contribution to the Ml decay width. However, because of the cancellation 
of the relativistic single quark current matrix element (^4|) for the cs system outlined in section 2.3, the 
two quark currents give the main contribution to the transition rate for D* s — > D s ^ p3[ . For the scalar 
confining interaction, that current may be expressed in the form 



j c {q,ki,k 2 ) 



Q*Pi Q2P2 



+ 0(^1+ ^2) X q 



-(oi -ct 2 ) X q 



Qi Q 2 



2m\ 

01 

2ml 



2m\ 



2ml 



(27) 



where q is the momentum of the emitted photon, which corresponds to p ' — p in case of the single quark 
amplitude. In eq. (p^), the variables Q\ and Q 2 are defined as Q\ = V c {k 2 )Q\ and Q 2 = V c (ki)Q 2 , where 
V c (k n ) is the (formal) Fourier transform of the scalar confining interaction. The momentum variables 
Pi and P 2 correspond to the expressions (pi +p' 1 )/2 and (p 2 +p' 2 )/2 respectively, while fc„ denotes 
the momentum transfered to quark n according to Fig. |l|. The analogous expression for the two-quark 
current induced by the vector coupled OGE interaction is of the form 



io\ x k 2 2P 2 + iS 2 x k 2 



2m\ 



2m\m 2 



ia 2 x k\ 2Pi + ia\ x k\ 



2m\ 



2m\m 2 



(28) 



with Ql = V g (k 2 )Qi and Q* 2 — V g (ki)Q 2 . Here V^(fc) denotes the momentum-space form of the OGE 
interaction. The corresponding magnetic moment operators may then be computed from eqs. (27) 
and (p8| ) according to 



/2=--lim V 9 x j (q,ki,k 2 ) 

2 q— >0 L 



(29) 



where the momentum transfer variables are given by k\ = q/2 — k and k 2 = q/2 + k. Upon Fourier 
transformation, the resulting magnetic moment operators may, for transitions between S'-wave states, be 
written in the form 



eV c (r) 



Qi Q2 



(<7l - <7 2 ) + 



Qi , Q2 



(ai + a 2 ) 



(30) 



for the scalar confining interaction, and 



eV g (r) 



Qi Q2 Qi — Q2 



m\m 2 



(a 1 - ct 2 ) + 



Q 2 , Qi 



m\m2 



{Si + <7 2 ) 



(31) 



for the one-gluon exchange interaction. In the limit mi = m%, the above magnetic moments reduce to 
those given in ref. ^7|. Note that for the case of equal constituent quark masses, the one-gluon exchange 
magnetic moment operator has no spin-flip term and consequently does not contribute to Ml decay if 
m\ = 777,2. The terms in eqs. ([30]) and ( |3l| ) that are symmetric in the quark spins affect only the magnetic 
moments of the heavy-light mesons and give no contribution to the Ml decay widths. 
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The spin-flip matrix elements for radiative Ml decay of a cs meson that should be added to the 
relativistic single quark matrix element are thus 



M 



Conf 

7 



dr u (r) 



V c (r) 



12 



1 



in case of the scalar confining interaction, and 



M Ogc = _ 



' dru 2 {r) Ml 
v ; 24 



m c m s 



(32) 



(33) 



for the one-gluon exchange interaction. In the above expressions, the appropriate values of the quark 
charge operators have been entered. For the scalar confining interaction, the form V c (r) — cr ~b will be 
used, where the parameters c and b are taken to be those obtained in the potential model of ref. |f23f . 
The one-gluon exchange potential V g (r) is taken to be of the form 



V g {r) 



4 2 
3?r 



dk a s (k 2 )j (kr) 



(34) 



which is known as the Richardson potential |B2j and features the running QCD coupling a s which is 
parameterized in terms of the QCD scale Aqcd and a dynamical gluon mass m g , which in ref. [p3| were 
obtained as 280 MeV and 240 MeV, respectively: 



a s (k 2 



12tt 



l 



27 ln[(fc 2 + 4m2)/A2 CD ] 



(35) 



Note that if a s is taken to be constant, then the above form reduces to the standard Coulombic one-gluon 
exchange potential of perturbative QCD. 



3.2 Two-quark contributions to ir° decay 

In ref. it was shown that a two-quark mechanism analogous to that for Ml decay shown in Fig. [j], 
which is due to coupling of the emitted pion to an intermediate negative energy quark, which interacts 
with the heavy spectator quark by the Qq interaction may contribute to the amplitude for pion decay. 
If the decay amplitude for emitted pions is written in the form 

T n = iq n ■ A, (36) 

then the axial current A may be decomposed into single- and two-quark contributions according to 
A = A s + A ex . The two-quark contributions, which have been calculated for equal constituent masses 
in ref. pq ], are in this case however proportional to m^ 3 in the nonrelativistic limit, which means that 
they are of much less importance than the two-quark contributions to the magnetic moment operator. 
Furthermore, it has been shown in ref. p6[ that the axial exchange current contribution associated with 
the scalar confining interaction is much reduced if relativistic effects are considered, because of the low 
mass of the strange constituent quark. 

The end result is, that two-quark effects due to the scalar confining and one-gluon exchange interac- 
tions will only modify the relativistic matrix element ( |2l| ) at the ~ 5 — 10% level. As the uncertainties 
in the value of the mixing angle 6 m and the empirically determined ratio of 7 and 7r° decay are much 
greater, there is no need to include the two-quark contributions to the amplitude for tt° decay at this 
time. 
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4 Numerical Results for D* s — ► £> s 7 and £>* — > D s ir° 

In order to obtain a description of the decay D* — » D s 7 that is as realistic as possible, the model for Ml 
decay presented here should be checked against the empirical results on the corresponding decays of the 
non-strange D mesons. As the decay width for D^* — > D ± j is known directly and that for D°* — ► D°7 
can be estimated from the empirical branching ratios and model calculations of the pionic decays of 
the D mesons, it is possible to calibrate the mass of the light constituent quark so that an optimal 
description of the Ml decay of the Z? ± is obtained. Together with the potential model calculation in 
ref. [^3| , this information can then be used to obtain a realistic estimate of the strange quark constituent 
mass. In order to calculate the Ml decay rate D ± * — > £ )± 7, it is sufficient to make the substitution 
m s — > m q in all matrix elements, but for D°* — > £>°7, one should substitute mj 1 — > —2m- 1 in eq. (p5|), 
and m" 2 — * —2m- 2 , l/m c m s — * 4/m c m q in eqs. ( |32] ) and (p3|). The obtained results for the Ml widths 
of the non-strange D mesons are given for different light constituent quark masses in Tables g and §. 



m q 


NRIA 


RIA 


RIA + Conf + Oge 


450 MeV 


0.58 keV 


9.4 


10- 3 


keV 


1.09 keV 


420 MeV 


0.79 keV 


1.5 


io- 2 


keV 


1.43 keV 


390 MeV 


1.07 keV 


2.2 


1Q- 2 


keV 


1.90 keV 



Table 1: Numerical results for the Ml transition D^* — > Z? ± 7, for a charm quark mass of 1580 MeV p3[ . 
The column NRIA (non-relativistic impulse approximation) corresponds to eq. (|2^), and the relativistic 
impulse approximation RIA to eq. (p4|). In the rightmost column, the two-quark contributions from 
eqs. (^2|) and ( [33] ) have been added. In the contribution from the confining interaction, the parameters 
c = 1120 MeV/fm and b = 320 MeV have been employed (23). 



m q 


NRIA 


RIA 


RIA + Conf + Oge 


450 MeV 


21.1 keV 


8.86 keV 


8.95 keV 


420 MeV 


23.5 keV 


9.18 keV 


9.89 keV 


390 MeV 


26.4 keV 


9.52 keV 


11.1 keV 



Table 2: Numerical results for the Ml transition D°* — > D°j, for a charm quark mass of 1580 MeV p3| . 
The columns NRIA (non-relativistic impulse approximation) and RIA are as for Table [lj although with 
the appropriate modifications due to the different quark charge operators. The same modifications apply 
to the rightmost column, where the two-quark contributions from eqs. ([50]) and (|l]) have been added. 
In the contribution from the confining interaction, the parameters c — 1120 MeV/fm and b = 320 MeV 
have been employed p3| . 
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It is instructive to compare the results presented in Table [y with the current experimental data for 
the D^* . Until recently, only the relative branching ratios for it and 7 decay were known p0[ . However, 
a first empirical measurement of the total width of the D^* has recently been published by the CLEO 
collaboration Q. The reported result is T(£) ± *) = 96 ± 4 ± 22 keV, where the latter error represents 
the systematic uncertainties. Taking into account the reported [^Oj branching ratio of 1.6 ± 0.4% for 
radiative decay, one obtains 1.5 ± 0.6 keV for the transition D^* — » D ±n f. Here most of the uncertainty 
stems from the systematic errors of the experimental result. However, it is evident that the results of 
Table [j] reproduce this result well for a range of values of the light constituent quark mass m q . The value 
ixiq = 450 MeV corresponds to the potential model of ref . (2j| , while the value m q — 420 MeV has been 
suggested by refs. |33[ |4|. From Table |l| it is thus seen that a light constituent quark mass of 420 MeV 
produces a width for radiative decay which is close to 1.5 keV. That value of the width is also favored 
by the analysis of ref. ]35| . 



m q 


NRIA 


RIA 


RIA + Conf + Oge 


560 MeV 


0.18 keV 


2.6 


10- 4 


keV 


0.38 keV 


530 MeV 


0.26 keV 


3.9 


10~ 5 


keV 


0.49 keV 


500 MeV 


0.36 keV 


8.6 


10- 4 


keV 


0.64 keV 



Table 3: Numerical results for the Ml transition Df* — > D^~/, for a charm quark mass of 1580 MeV p3[ . 
The column NRIA (non-relativistic impulse approximation) corresponds to eq. (|2^) , and the relativistic 
impulse approximation RIA to eq. (|24j). In the rightmost column, the two-quark contributions from 
eqs. ([[2]) and (|33] ) have been added. In the contribution from the confining interaction, the parameters 
c = 1120 MeV/fm and b = 260 MeV have been employed [^3|. Note that the RIA contribution becomes 
exactly zero for a strange quark mass of ~ 540 MeV. 

Even though the total width of the neutral D°* meson has not yet been determined considerable 
information about the expected width for D°* —> D°j may be extracted from the reported branching 
fraction of 38.1 ± 2.9% |^0| , since the corresponding width for pion decay can be constrained by means of 
the empirically determined width of the Z? ± * and model calculations of the pionic decays of D mesons |i~6[ 
|l7j , f36[ . This is seen directly if one first notes that the branching fraction of Z? ± * — >• Z?°7r ± is reported as 
67. 7± 0.5% |2Cj , which implies a width for this decay mode of ~ 65 ± 14 keV. From the model calculation 
of ref. |l(| one finds that this corresponds to a width of ~ 40 ± 10 keV for D a * — > D°ir°. As the relative 
branching fractions for tt° and 7 decay of the D°* are well known J2(J , the best estimate for the width 
of D°* — > Z?°7 is ~ 25 keV. There remains, however, still a considerable uncertainty of ~ ±10 keV from 
the systematic errors in the empirical measurement of T(D ± *). 

A width for D°* — > of around 20-30 keV is also preferred by ref. (35| . The result of the model 
calculation in Table || is smaller by about a factor ~ 2 as the OGE contribution is weaker in this case as 
compared to the contribution from the scalar confining interaction. Moreover, as a relativistic treatment 
of the two-quark currents will effectively lead to replacement of the quark masses by energy-dependent 
factors, it may be argued that relativistic effects will generally lead to a weakening of the two-quark 
contributions to the matrix element for Ml decay. It is thus entirely possible that a fully relativistic 
treatment will render the two-quark contributions too weak to account for the experimental data on Ml 
decay of the Z? ± meson as well. 
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This conclusion is in line with that reached in ref. [ pq ], where it was found that the introduction 
of a large anomalous magnetic moment is required in order to fit the observed Ml decay widths of 
the heavy-light mesons. Another interesting possibility is that the instanton induced interaction for Qq 
systems given in refs. |3^, |34j may also contribute a significant two-quark current. The interaction of 
refs. (3^, [54|, which was found to be short-ranged, negative and attractive, has scalar coupling to the light 
constituent quark. One may thus conclude that this interaction will add up constructively with the OGE 
contribution from eq. (33). Thus the instanton induced interaction will have an overall favorable effect 
on the widths for Ml decay, particularly for the D *, which may be inferred from the matrix elements 
given in Table f|. 

As it appears to be possible to obtain a realistic description of the radiative Ml decay of the Z? ± 
meson using the model presented in this paper, then a prediction based on the potential model p3| for 
the Ml decay Df* — > Dfj and the associated matrix element A4 y can be made. That information is 
needed as input for eq. (|13|) in order to obtain an estimate for f nss - These results are given in Tables | 
and |. 





D ±* D ± 1 


D°* -> D a 1 




M Rcl 

M° oni 
M° sc 


-1.83- 10~ 2 
-2.55 • 1(T 3 
1.23- 10" 2 
-3.45 • 1(T 2 


9.91 • 1(T 2 
6.20- 1(T 2 
-3.07- 1(T 2 
3.31 • 1(T 2 


-8.55- 1CT 3 
3.24- 1CT 4 
7.25 • 10^3 

-1.98- 1CT 2 


M Tot 


-2.48- 1(T 2 


6.44- 1(T 2 


-i.22 • itr 2 



Table 4: Matrix elements for the Ml decays considered in Tables [j], ||and |[ in units of [fro]. The matrix 
elements correspond to a charm quark mass of 1580 MeV, a light constituent quark mass of 420 MeV 
and a strange constituent quark mass of 560 MeV. The matrix element A4^ ot shows the sum of the 'Rel' 
+ 'Conf ' + 'Oge' contributions, which in case of the decay Df* — ► Dfj is used for the determination of 
the coupling constant f V ss- 



In order to estimate the coupling constant , the matrix element (^) for n° decay of the D* meson 
also needs to be evaluated. The value of that matrix element for a strange constituent quark mass of 
560 MeV is 0.794, which together with the total matrix element for 7 decay of the D* may then be used 
with eq. (^) to obtain a value for f^ss- If the charm quark also couples significantly to the 77 meson, 
then the matrix element for pion emission by the charm quark in eq. (15) needs to be evaluated as well. 
This can be done by the substitution m s — * m c in eq. (|2l]). For a charm quark mass of 1580 MeV, the 
result M.% — 0.949 is obtained. Thus, in the evaluation of eqs. ( |l4| ) and (15), the following values for 
the matrix elements for 7 and ir° decay will be used: 



= -1.22- 10" 2 fm 



for the 7 decay matrix element, and 
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M% = 0.794 M% = 0.949 



for the 7T° decay matrix elements. Note that in all the calculations of this section, the masses of the 
initial and final quarkonium states as well as the momenta of the emitted photons and ir° mesons have 
been taken to equal those given in ref . |Q ■ 



4.1 Numerical Results for ip' — > J/ipr] 

The evaluation of eq. ([l6j) involves the computation of the matrix element A4 V given by eq. (^2|) . The ip' 
and J I ip wavefunctions are here taken to be those computed in the potential model of ref. fl22[ . Because 
of the somewhat larger value of q v (200 MeV) pPf , and the orthogonality of the ip' and J/ip radial 
wavefunctions, it is undesirable to make the approximation q n = in this case. The resulting relativistic 
matrix element for the decay ip' — > J /ip rj is obtained as A4 V = —3.53 • 10~ 2 , and gives the desired 
width of ~ 7.5 keV if the coupling constant f vcc is taken to be ~ 0.8. The sign of f vcc is of course not 
determined by this computation. 

As there exists both theoretical and phenomenological arguments js] against a significant cc admixture 
in the rj meson, then this value for f^cc may appear uncomfortably large. Indeed, it will be shown in 
the next section that a value of f vcc around ~ 0.8 is, in view of the D* — ► D s tt° decay, about one half 
of the so obtained value for f vss - However, as neither the axial charge component nor the associated 
two-quark contributions to the amplitude for r\ decay of the ip' have been considered, the value of 0.8 
for f vcc should be considered as an upper limit only. Nevertheless, it will be shown in the next section 
that the inclusion of a pion-charm quark coupling will serve to increase the obtained value for f vss , thus 
bringing it closer to the phenomenologically preferred values 
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5 Discussion 



Assuming that in the ir° decay of the D*, the pion couples mostly to the strange constituent quark, 
eq. (|l3| ) may be used directly together with the matrix elements for ir° and 7 decay given in the previous 
section to estimate the magnitude of f V ss- Insertion of the matrix elements and taking 9 m = 0.012 yields 
R = 0.3085, giving |/, ;ss | = 0.70. Taking into account the uncertainties in the mixing angle and the 
empirical decay widths for 7r° and 7 decay, the best estimate is 



fr]s 



0.7 



+0.5 
-0.3 



In the above result, the negative sign is suggested by the relations in eq. (g). The static quark model 
then implies, through eq. (Q), that the magnitude of the corresponding pseudovector 77-nucleon coupling 
constant f V NN should be one half of this value. Thus one obtains the following value for the ry-nucleon 
coupling: 



jjNN 



= 0.35 



+0.15 
-0.25 



This result should be compared with the value for f V NN (or the equivalent pseudoscalar coupling constant 
9t)NN = (2mjy /m n )f V NN , which has been determined by phenomenological model fits to photoproduction 
of 77 mesons on the nucleon jt^]. The latter value for f V NN is 0.64. This value has also been found to 
be realistic in calculations of the cross section for pp — > ppij near threshold Jll|. Although the result 
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obtained above for fr/NN has quite large uncertainties which are mostly of empirical origin, it still appears 
to be significantly smaller. A larger value for f^NN could of course be obtained by decreasing the mixing 
angle 9 m . 

Additional modifications to the current estimation of f n NN may of course also arise from SU(3) 
breaking effects in eq. (H), which relates f mq to f vss . Likewise, relativistic effects may also modify the 
static quark model result (Q) which relates f mq to f V NN- However, in view of the large experimental 
uncertainty in the determination of the branching fractions for 7 and 7r° decay of the D* , there appears 
to be little motivation to include such effects at this time. Another possibility is that the matrix element 
for 7 decay as given by Table f| is too small. However, as the computed width for D* — > D s ^ is already 
somewhat larger than that given in ref. |3j|, this possibility may be considered unlikely. 

Another option worth considering, in view of the fact that the empirical detection of a considerable 
branching fraction for tj) 1 — > J ftp 77 indicates that the 77 meson also couples to the charm quark, is that 
the coupling of the 77 meson to the charm quark would be strong enough to significantly influence the 
D* —>■ D s ~f decay. For nonzero values of the coupling f-qcc, the coupling constant f vss may thus be 
determined from cq. (|l5| ) according to 

_ VR- f vcc M% 

hss- M% , 

where the different possibilities are y/R = ±0.5554 and f^cc — ±0.8. The most realistic of these appears 
to be \/R = —0.5554, f vcc = +0.8, giving f vss = —1.66. In this case a nonzero and positive /^ cc leads 
to a larger and negative value for f VS s> allowing for better agreement with the values of f n NN suggested 
above. The value for fr/NN that corresponds to f ncc — ±0.8 turns out to be 

f V NN = 0.83 (38) 

which is closer to the phenomenological estimates given by refs. |Q and Although the value ±0.8 
for frjee is, in view of the discussion in the previous section, probably much too large, there still remains 
a distinct possibility that the charm quark contributes significantly to the width for D* — > D s ir°. In this 
case it is found that values around /^ ss ~ —1.7 are favored by an eta-charm coupling of ~ +0.8. 

It is useful to compare the results obtained from the present phenomenological analysis with the 
predictions of eq. (^). With the value g q A — 0.87 for the axial coupling constant of the constituent 
quarks, eq. (Q) yields the value f^ qq — 0.65 for the pion-quark coupling, which is close to the value 
f-nqq = S/bf^NN — 0.57 suggested by the static quark model (cf. @). The values for f vqq and /, )ss so 
obtained are 1.25 and -2.5 respectively. The values for f vss obtained here from the ratio of the empirical 
branching fractions for D* — * D s ir° and D* — * D s j, are much smaller than this quark model value. In 
the completely SU(3) flavor symmetric limit, in which the pion and rj are degenerate in mass, the value 
for f, lss would however be equal to —2/v3fnqqi which is only -0.75. The results from various modern 
phenomenological analyses are typically about twice this value. 

Naive application of SU(3) flavor symmetry at the level of a chiral symmetry violating pseudoscalar 
coupling model for pions and r\ mesons to baryons gives the result g V NN = 1 / g^NN- With the (now) 
standard value of 12.8 for g^NN this relation would give the value 7.4 for g V NN- Such a large value would 
correspond to J v nn — 2.2. While values of that magnitude are ruled out by 77 meson photoproduction, 
as well as, as shown here, by the 7r° decay of the D* meson, they were employed in early realistic boson 
exchange models for the nucleon-nucleon interaction. 
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